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Abstract 



We discuss the nonmesonic decay of deformed p-shell hypernuclei. The 
Nilsson model with angular momentum projection is employed in order to 
take into account the deformation effects. The nonmesonic decay rate and 
the intrinsic A asymmetry parameter decrease as a function of the deforma- 
tion parameter, while the ratio of the neutron- to proton- induced decay rates 
increases. We find that the deformation effects change these observables by 
about 10 % for ^Be from the spherical limit. 



1 



Typeset using REVT^ 



I. INTRODUCTION 



One of the main issues of nuclear physics is to understand the nucleon-nucleon (NN) 
interaction. The \AS\ = 1 NN interaction is particularly important in this respect, since the 
change of strangeness can be used as a signature to study both the parity-conserving (PC) 
and the parity-violating (PV) amplitudes. This is in clear contrast to the AS = weak NN 
interaction, where the weak PC signal is masked by the strong interaction. 

Due to the lack of stable A-particle beams, the weak decay of A-hypernuclei has been 
the only source of information on the weak four-baryon \AS\ = 1 interaction. Single A- 
hypernuclei are typically produced via either hadronic reactions, as {Kg^gp,7i'^) ^ or {7i^,K^) 
0, or electroproduction mechanisms, as (e,e' K^) These hypernuclei are typically 
produced in an excited state and reach their ground state by electromagnetic-7 and/or 
particle emission. Once they are stable against strong decay, they decay via weak interaction 
mechanisms which are nonleptonic in nature and violate isospin, parity and strangeness. 
Since the mesonic decay mode, A ttN, is Pauli blocked in the nuclear medium, hypernuclei 
with A> 5 predominantly decay through the nonmesonic decay (NMD) mode, AA^ NN. 

In order to learn about the weak AA^ NN interaction from the theoretical side, one 
has to take into account different inputs as accurately as possible. These include the de- 
scription of nuclear structure, the choice of the strong BB potential model [^,0, A/ =1/2 
violations P] and the importance of the 3N emission channel, Anp nnp In Refs. 

PJTOd, a one-meson-exchange (OME) model was applied to calculate the nonmesonic decay 
observables of the p-shell and and the s-shell |He and hypernuclei. We included 
the virtual exchange of the ground-state pseudoscalar and vector mesons p, ri,u, K and fT*, 
in addition to the long-ranged pion. Except for the hypertriton, where the hypernuclear 
wave function was calculated exactly using the Faddeev formalism, the structure of the 
initial hyper nucleus was described in a shell-model framework which assumed spherical con- 
figuration. In these calculations, the strong baryon-baryon (BB) interaction was accounted 
for using the Nijmegen BB potential model H]. Monopole form factors at each vertex were 
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included in order to regularize the weak potential, while the weak baryon-baryon-meson 
coupling constants were derived based on SU^(6) and soft-meson theorems. The total NMD 
rate and the asymmetry in the distribution of emitted protons from the decay of polarized 
hypernuclei were in good agreement with the experimental data. However, the theoretical 
values for the neutron-to-proton ratio were found to be very small compared to the exper- 



imental data. Several attempts have been made to reconcile this discrepancy [pHgJTlHlS 
but none of them has solved this problem yet. 

Our aim in this paper is to investigate how much these observables depend on the de- 
formation of hypernuclei. All previous calculations were performed using the spherical con- 
figuration, however, it is well known that many p-shell nuclei are deformed in the ground 
state. For instance, the quadrupole deformation parameter extracted from the experimental 
quadrupole moment [1^ is (32 = 0.65 for and —0.71 for ^^C. It may be important to take 
these deformation effects into account in order to describe quantitatively the nonmesonic de- 
cay of p-shell hypernuclei. Deformed hypernuclei can be described using several models such 
as the a-cluster model []T3[ or the deformed self-consistent Hartree-Fock method. In fact, one 
can also use realistic wave functions obtained by a diagonalization of a shell-model Hamil- 
tonian for p-shell nuclei, as in Ref. []TB[. In the present paper, however, in order to perform 



a systematic study, we use instead the Nilsson model |117ip!8| as a simplified Hartree-Fock 
method. 

The paper is organized as follows. In Sec. II, we present the relevant formulae to 
evaluate the NMD observables in a OME model. In Sec. HI, we briefly review the deformed 
shell model based on the Nilsson model. Sec. IV presents the deformation dependence of 
the nonmesonic observables for the decay of \Be, whose ®Be core is known to be largely 
deformed. Although there is no experimental data for this hypernucleus at present, we 
choose this system as the simplest non-spherical p-shell hypernucleus and as a representative 
example of deformed p-shell nuclei. We compare our theoretical predictions with the typical 
experimental data for other p-shell hypernuclei. Sec. V summarizes the paper. 
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II. NONMESONIC WEAK DECAY IN A ONE-MESON-EXCHANGE MODEL 



Assuming that the initial hypernucleus is at rest, the NMD rate is 

{1}{2} 

where A^/i is the hypernuclear transition amphtude. The quantities Mh, E^, Ei and E2 
are the mass of the hypernucleus, the energy of the residual {A — 2)-particle system, and the 
total asymptotic energies of the emitted nucleons, respectively. The integration variables ki 
and k2 are the momenta of the two baryons in the final state. The momentum conserving 
delta function has been used to integrate over the momentum of the residual nucleus. The 
sum, together with the factor 1/(2 J + 1), indicates an average over the initial hypernucleus 
spin projections, M/, and a sum over all quantum numbers of the residual {A — 2)-particle 
system, {R}, as well as the spin and isospin projections of the exiting particles, {1} and {2}. 
In general, one can write the total nonmesonic decay rate as = r^^^^^ = Fn + Fp, 
where Fn {An nn) stands for the neutron-induced decay and Fp {Ap np) for the proton- 
induced one. 

In addition to the total and partial decay rates, we also calculate the intrinsic A asym- 
metry parameter. When working with polarized hypernuclei and in combination with co- 
incidence measurements of the decay particles, one can study the angular distribution of 
particles coming from the AA^ NN weak decay. Due to the interference between the PV 
and PC amplitudes, the distribution of the emitted protons in the weak decay displays an 
angular asymmetry with respect to the polarization axis. The asymmetry A, defined by 

_ 3 TrjMf.SyMf.^) 

is expressed in terms of the hypernuclear polarization created in the strong production re- 
action, Py, the J-spin operator along the polarization axis, 5*^^, and the total spin of the 



initial hypernucleus, J. In Ref. [O] it is shown that the asymmetry follows a simple cosx 



dependence, i.e., A = PyAp cos x, where x stands for the angle between the direction of the 



proton and the polarization axis. The hypernuclear asymmetry parameter Ap is character- 
istic of the hypernuclear weak decay process and depends on J and the intensity of protons 
exiting along the quantization axis for the different spin projections of the hypernucleus. 
At X = 0°, the asymmetry in the distribution of protons is thus determined by the product 
A = PyAp. In the following, we assume a weak coupling scheme where the A hyperon is 
coupled only to the ground state of the {A — l)-particle core. In this scheme, simple angular 
momentum algebra relates the hypernuclear polarization Py to the A polarization p\, and 
the hypernuclear asymmetry parameter Ap to the intrinsic A asymetry parameter a a, such 
that A = PaOa = PyAp . 

The nonmesonic decay of hypernuclei proceeds through a two-body mechanism. There- 
fore in order to evaluate the transition amplitude in Eq. (|l]), one has to decompose the 
(A-l)-core wave function into a set of states in which a nucleon couples to the residual (A- 
2)-particle state. This can be done using the Coefficients of Fractional Parentage (CFP), 
which are defined by 



where and Tr are the spin and isospin of the residual nucleus. The weak potential re- 
sponsible for this transition can be obtained by making a nonrelativistic reduction of the 
free Feynman amplitude depicted in Fig. [1|. In Table | we show the strong and weak vertices 
for pseudoscalar (PS) and vector (V) mesons. A,B,a,P and e stand for the appropriate 
baryon-baryon-meson weak coupling constants, while g {g^ , g"^) represents the strong (vec- 
tor, tensor) coupling. Details of the derivation of the transition potential can be found in 
Ref. and here only the final expression will be presented. For pseudoscalar mesons, the 
potential is 



where Gf^-k"^ = 2.21 x 10 ''is the Fermi coupling constant, g is the momentum carried by 
the meson directed towards the strong vertex, /i the meson mass and M (M) is the average 



JM,TT,)= J2 {JT{\jRTn,jt) | J^jT^) (g) | jt) 



J JM,TT^ 



(3) 
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of the baryon masses at the strong (weak) vertex (the other way around for the exchange of 
strange mesons). For vector mesons, the potential is 



K(g) = Gprn^ \ g^a (cti x q){a2 x q) 

'(ai xa2)g ^— . (5) 



The values of the strong and weak couplings are listed in Table III of Ref. [^. In Eqs. (^ 

and (^ the operators A, B, a, (3 and e contain, apart from the weak coupling constants, the 

specific isospin dependence of the potential, which is fi ■ T2 for the isovector vr and p mesons, 

i for the isoscalar rj and uj mesons and a combination of both operators for the isodoublet 

K and K*. In order to derive Eqs. (^) and (^ we assumed the validity of the A/ = 1/2 

rule, which is known to experimentally dominate the decay of A's into pions. A/ = 3/2 

transitions for vector mesons (p and K*) are easily accomodated P] in the formalism and 

the results we present here account for such A/ =1/2 violations. 

We obtain a regularized potential by including a monopole form factor at each vertex, 
(A^ - p2) 

Fiq"^) = TTT, where the value of the cut-off. A, different for each meson, is taken 

(A2 + g2)' 

from the Jiilich hyperon-nucleon interaction [Q . To incorporate the effects of the strong NN 
interaction, we solve a T-matrix scattering equation in momentum space for the outgoing 
nucleons using the Nijmegen potential models. For the initial bound two-baryon system 
we use a spin independent parametrization of the type 

/(r) = fl - e-"'/'^')" + for^e-'-'/^' , (6) 



with a = 0.5 fm, b = 0.25 fm~^, c = 1.28 fm, n = 2. The results obtained with this 
parametrization pO[ lay in between of the ones obtained with a microscopic finite-nucleus 
G-matrix calculation [pl|] using the soft-core and hard-core Nijmegen |23] models. 



III. DEFORMED SHELL MODEL 



As we mentioned in the previous section, we use a weak coupling scheme for the A hyperon 
in the initial hypernucleus. To this end, we must describe the ground state of the core nucleus 
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which may be deformed. The Nilsson model provides a simple and convenient framework 



to describe deformed nuclei, and has been widely used in the literature [p^ , p^ , p^ -p6[] . Its 
Hamiltonian consists of an anisotropic harmonic oscillator with the spin-orbit interaction as 
well as an angular momentum dependent term, which mocks up the deviation of the mean 
field potential from the harmonic oscillator 

H = Ho-^^6muj',Y,oie), (7) 

= -^V' + \mulr'' + Cl- s + D {f - {f)^) - 6 mujlY^oie). (8) 

Here 5 is a deformation parameter, / and s are the single-particle orbital and the spin angular 
momenta, and C and D are adjustable parameters. (/ )iy = N{N + 3)/2 is the expectation 
value of / averaged over one major shell with quantum number N. The relation between 5 
and P2 is given by 



Since the Nilsson Hamiltonian (H) violates rotational invariance, the total angular momen- 
tum j = I + s is not a good quantum number. However, the projection of j on to the 
2;-direction, k, is conserved and the single-particle levels are characterized by k and other 
quantum numbers. We expand a Nilsson single-particle level, ipk{q), in terms of the eigen- 
functions of the spherical harmonic oscillator Hamiltonian Hq, (pnijk, as 

= J2^nljk4>nljk, (10) 
nlj 

where q are quantum numbers other than k. We choose x^^ij_i^ = {—y~''x^nijk ^'^ ^^^^ 
eigenvalues of the Nilsson Hamiltonian do not depend on the sign of the projection of total 
angular momentum We denote the creation operator of ipk as a|. and that of (pjk as 6^^. 



We explicity express only the j and k quantum numbers to simplify the notation. Intrinsic 
wave functions, i.e., eigenfunctions of the Nilsson Hamiltonian (H) are given by 

\Xk) = alal ■ ■ -aljO) = (y^x^k^k) |0), (H) 



where the K quantum number is the sum over all ki. The intrinsic wave function ([TT|) is not 
an eigenstate of the total angular momentum J, and thus has to be projected out to a good 



angular momentum state. This can be achieved by using the projector given by p8| , p3| , p4 

2J+ 1 



MK 



in" 



(12) 



where Q are Euler angles, and Dfjj^lQ) and R{Q) are the Wigner D-function and the rotation 
operator, respectively. 

For systems with a single-Nilsson level, such as ^Be which we discuss in the next section, 
the CFP can be analitically obtained [^,^ . Note that a single Nilsson level can accommo- 
date up to four nucleons, i.e., two protons and two neutrons. For three particle systems ( 2 
neutrons and 1 proton, for example), the wave function is given by 



^JM = [N{3)j]-'Pi,^j,^,aUUAjO), 

= [A^(3)j]"^ ^ y^^Xj^kXj2-kXj3k{jikj2 - k\Ji20){Ji20j3k\Jk) 

ilJ2,i3 Jl2 



(13) 



X 



|0), 



(14) 



JM 



where tt stands for proton and z/ for neutron. The normalization factor N{?>)j is given by ||23 



[A^(3)j]' = $^5j,,,e.en(x,fc)'f/(Ji2,A;)(Ji20jA;|JA;)^ 

3,Jl2 



(15) 



where 



f/(J,fc) = 2Y,{^nk?{^32k?{3ik32 - k\.m). 
The isospin of this system is 1/2. For four nucleon systems, the wave function reads 



(16) 



(17) 



[N{A)j] ^ X] X] Xj^kXj.,-kXj^kXj^-k{iiki2 - k\Ji20){Ji20j3k\Jiik){Jiikj4 - k\JO) 

jl,j2 h,3i Ji2,Jr 



«k 1 |o)' 

•^R / JM 



with the normahzation given by 

[^(4)j]'= Sji2,even5j,„evenU{Jl2,k)U{J3^,k){Ji20J340\J0)\ (19) 

The isospin of this wave function is 0. Comparing Eqs. (|I4D and (|18D, the CFP for the four 
particle system reads 

C{j) = {JT{\JnTn,jt) = -v^x,_,(J«fcj - A;| JO)^^^. (20) 

[4:) J 



IV. NONMESONIC DECAY OF ^BE 

Let us now apply the deformed shell model of Sec. Ill to the nonmesonic decay of \Be. 
The quadrupole moment of the neighbour nucleus ^Be was measured to be 5.86 efm^ |]14 



from which we extract the quadrupole deformation parameter /32=1.00 using the radius 
parameter ro=1.2 fm. Several theoretical calculations suggest that the core nucleus ®Be and 
the ^Be hypernucleus also have similar deformation parameters with the same sign |]T5[. Our 
interest is to discuss such deformation effects on nonmesonic decay observables. 

As is discussed in Sec. II, the use of CFP allows us to write the hypernuclear transition 
amplitude Aifi in terms of elementary two-body amplitudes. Therefore, our first task is 
to compute these coefficients for the core nucleus '^Be. We assume the inert spherical "^He 
core and explicitly work with only the four valence nucleons. Diagonalizing the Nilsson 
Hamiltonian (|^) , one finds that the lowest Nilsson level for the valence nucleons has k = 1/2 



for prolate deformation [Q. We diagonalize the Nilsson Hamiltonian in the AA^=0 states. 
Contributions from the AA^ = 2 can be neglected unless the deformation is large. The 
k=l/2 state is thus 

where x and y are determined by diagonalizing the Nilsson Hamiltonian within this configu- 
ration space and depend upon the deformation of ^Be. Using Eq. (pUD, the CFP's are found 
to be 



\C( = 3x« + ?,xW + .22^ 

for the P3/2 state, and 

^^2 _ x^y"^ + 9a;^2/^ + 

for the pi/2 state. Note that [C(]93/2)]^ + [C(pi/2)]^ = 1. In the spherical hmit, x=l and 
7/ = 0, so the CFP become [C(]93/2)]^ = 1 and [C(pi/2)]^ = 0. The CFP for the deeply bound 
lsi/2 state is just equal to 1 since ^He is a spin-isospin saturated nucleus. 

Our results for the nonmesonic decay rate, F"™, in units of the free A decay rate. Fa = 
3.8 X 10^s~^, the neutron-to-proton ratio, Fn/Fp, and the A asymmetry parameter, oa, are 
shown in Fig. 2 as a function of the deformation parameter (32- We use an oscillator length 

of 1.65 fm for nucleons, so that the experimental root mean square radius of ^Be is 
reproduced. Following Refs. [0,IT8[, the parameters C and D in the Nilsson Hamiltonian 



(1) are taken to be — O.lG/kuo and 0, respectively. As for the oscillator length h\ for the 
lsi/2 wave function of the A hyperon, we estimate it to be 1.5 fm in order to reproduce its 
binding energy in ^Be (= 6.71 ± 0.04 MeV [0). From the figure, we see that F°™/Fa is a 



decreasing function of [32, while Fn/Fp and oa are increasing functions. As we have already 
mentioned, the deformation parameter of ^Be is expected to be close to 1. We notice that the 
nonmesonic decay observables are altered by about 10% from the spherical limit at /52 = 1- 
An important question is whether this effect is significant when comparing to the exper- 
imental data. We note that the typical experimental uncertainties for nonmesonic decay of 
p-shell hypernuclei are: 7% - 17.5% for the total decay rate [g|Hli' 46.2% - 84.2% for the 



neutron-to-proton ratio [p9| , p0| , and 50% - 1000 % for the asymmetry [^. These experimen- 
tal uncertainties are much larger than the theoretical one originating from the deformation 
effects. Thus we conclude that the spherical approximation gives a good estimate of the 
nonmesonic decay of p-shell nuclei, at least within the present experimental precision. 

Before we close this section, we would like to stress that our conclusion is not altered 
qualitatively even if more realistic wave functions are used instead of the present schematic 
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ones. For instance, using the shell-model C.F.P. of Cohen and Kurath |T6| for the decay of 



]^C, we obtain F^^/Fa = 0.76 and Fn/Fp = 0.22. Those numbers have to be compared to 
the spherical limit values of 0.74 and 0.25, respectively. As we see, the amount of deviation 
of those observables with respect to the spherical limit is of the same order of the one given 
here, although their behaviour is opposite due partly to the fact that the residual ^^C is 
oblate while *Be is prolate. 

V. SUMMARY 

We have discussed the role of nuclear structure in the nonmesonic decay of p-shell hy- 
pernuclei, especially focusing on the effects of deformation. To this end, we have used the 
Nilsson model with explicit angular momentum projection. We have studied the nonmesonic 
decay of ^Be as a typical example of deformed p-shell hypernuclei. We have shown that 
the deformation effects change the total NMD rate, the neutron-to-proton ratio and the A 
asymmetry parameter by about 10 % from the spherical limit. Although this value is not 
negligible, it still is smaller than the present typical experimental uncertainty and smaller 
than other theoretical uncertainties, e.g., the effects of SU(3) symmetry breaking ||9|,p!2 



or 



AI =1/2 violations [Q. This indicates that the existing discrepancies between the experi- 
mental and theoretical values of hypernuclear weak decay observables cannot be attributed 
solely to deviations from the spherical configuration, and still remain an open question. New 
experiments are urged in order to reduce the large experimental error bars, which prevent 
any definite conclusion about the reliability of the theoretical models. 

Our conclusions may not be the same for heavier hypernuclei such as ^^^U ||32| , |33| . There 
are a lot of intruder states in such heavy deformed systems, unlike p-shell nuclei where there 
is only a few, or maybe zero, intruder states. Therefore, an interesting future work would 
be to discuss the nonmesonic decay of heavy hypernuclei including the deformation effects. 
For that purpose, the projected shell model developed in Refs. p3| , p^ , which also uses the 
Nilsson model with angular momentum projection, would provide a powerful tool to describe 
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the structure of deformed hypernuclei. 
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FIGURES 

FIG. 1. Free Feynman diagrams for the AN — ^ NN transition mediated by the exchange of 

the nonstrange vr, rj, p, to (left) and strange K, K* (right) mesons. The shaded circle (filled) stands 
for the weak (strong) vertex. 
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TABLES 

TABLE L Vertices entering the expression of the Feynman amphtude of Fig. |^ for pseudoscalar 
(PS) and vector (V) mesons. The weak vertices are in units of Gpiri.^'^. 
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